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We study the Berezinskii-Kosterlitz-Thouless transition of two-dimensional Bose gases in a syn- 
thetic magnetic field using the standard Metropolis Monte Carlo method. The system is described 
by the frustrated XY model and the critical temperature is calculated though the absence of central 
peak of the wave function in momentum space, which can be directly measured by the time-of- flight 
absorbing imaging in cold atoms experiments. The results of our work show agreement with former 
studies on superconducting Josephson arrays. 
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I. INTRODUCTION 

It is well known that in two-dimensional (2D) sys- 
tems with a continuous symmetry the conventional long- 
range order is prevented by thermal fluctuations at fi- 
nite temperature in the thermodynamic limit, and as 
a result no spontaneous continuous symmetry breaking 
takes place. However, these systems can undergo a tran- 
sition through binding of vortex pairs to form a quasi- 
long-range order. This is the celebrated Berezinskii- 
Kosterlitz-Thouless (BKT) transition 0, 0. Theoreti- 
cally, the 2D XY model is the prototype to elaborate on 
the BKT transition. Experimentally, the BKT transition 
has been examined in various physicsl systems, including 
4 He films Q, 2D superconductors Q and superconduct- 
ing Josephson-junction arrays (JJA) [Ej], and experimen- 
tal works agree well with theoretical predictions. 

The study of Josephson-junction arrays greatly en- 
riched the research of the BKT transition [5|42[- Once 
a transverse magnetic field is introduced to JJA, the 
BKT transition point as a function of the field can be 
investigated [5, 6]. Theoretically, such a system can be 
described by the frustrated 2D XY model 0-0 • Never- 
theless, these systems reveal a complex energy spectrum, 
known as the "Hofstadter butterfly" [Toj| , which gives rise 
to rich phenomena that are still of present research in- 
terest. 

So far, cold atoms have been regarded as ideal test beds 
for fundamental models of condensed matter physics. It 
is of great interest to examine the BKT transition since 
the quasi-2D Bose gas has already been produced ei- 
ther in a single pancake traps or at the nodes of one- 
dimensional optical lattice potentials Trombettoni 
it et al. has raised a proposal to investigate the occur- 
rence of the BKT transition in a 2D scalar Bose-Einstein 
condensate [12|. The BKT transition in atomic Bose 
gases was first experimentally realized in 2006 by Hadz- 
ibabic [l3j|. They added the optical lattice in the z di- 
rection to construct a quasi-two-dimensional configura- 
tion, and the transition was detected by the dislocations 
of interference patterns. Later experiments identify this 
transition by the proliferation of vortex pairs [14| and a 



bimodal density distribution or coherent length [~L5^ . 

The newly developed technic of synthetic magnetic 
field in cold atoms use the light-modulated excitation 
between nearby neighbor sites to construct a "magnetic 
field" like effective hamiltonian in a 2D optical lattice 
[l6|. This technic opens the door for cold atoms to sim- 
ulate some more profound models, e.g. the lattice gauge 
theory. Thus it stimulates us to re-examine the behavior 
of BKT transition in the cold atomic Bose gases subject 
to the synthetic magnetic field. Comparing with the for- 
mer 2D experimental systems concerning on the super- 
fluid density or sheet resistance to detect the superfluid 
transition, cold atoms can be manipulated more easier 
and precisely and thus it benefits for further investiga- 
tion of the BKT transition. 

We start with the Bose-Hubbard model which de- 
scribes a 2D Bose gas in uniform magnetic field in a 
square optical lattice. This model can be mapped to the 
frustrated XY model. Then we study the BKT transi- 
tion using the standard Metropolis Monte Carlo method. 
This method has been used by Trombettoni et al. to 
investigate the BKT transition in a atomic Bose gases 
in a 2D optical lattice without magnetic field [12]. The 
transition point is determined through detecting the cen- 
tral peak of the expansion condensate after time-of-flight 
(TOF) expansion. 



II. THE MODEL 

Consider a 2D Bose gas immersed into a uniform mag- 
netic field in a square optical lattice. The Hamiltonian 
of this system can be described by the frustrated Bose- 
Hubbard model: 

fi = - K Y, (^Ujz iAij + h - c ) + ^ - 1) , (!) 

(id) 

where the first term on the right side is the hopping en- 
ergy with the summation over the nearest neighbor sides 
and K is the hopping matrix element. ^(^J) is the field 
operator of bosons for annihilate (create) a boson at site 
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i. Aij is a bond operator describing the magnetic field 
and around every plaquette we have: 

Aab + A bc + A cd + A de = 2tt/. (2) 

Here / = j^Hoa 2 = ^ is the uniform frustration with 
0o = ftc/2e the flux quantum and we have chosen unit 
lattice constant a = 1. The second term of Eq. (pQ) refers 
to the on-site repulsive (attractive) interaction, according 
to the interacting strength U > 0(U < 0), and Ni is the 
particle number operator on site i. 

The system under consideration is a uniform system at 
sufficient low temperature, such that the system is den- 
sity coherent with fluctuating phases. The field operator 
of bosons can be written as ipi = y^Noe %ei , where Nq is 
the average particle number of each site and 0i is the 
corresponding phase on site i. The mechanism of BKT 
transition, the pairing of vortices, is due to topological 
long-range correlation, thus we can safely ignore the on- 
site interaction term at sufficient low temperature. Upon 
these, the Eq. (pp) can be mapped to the frustrated XY 
model: 

H = -JY J cos ( i ~ °J + A v) ■ ( 3 ) 

(iJ) 

Here J = 2KNq. This model was first used by Titel et 
al. to describe the superconducting Josephson arrays in 
transverse magnetic field Q. 

We note here that the frustrated XY model is U(l) 
gauge symmetry breaking in the superfluid state as the 
hopping term breaks the conservation of " charge" . It is 
to say that in a conventional system described by the 
frustrated XY model physical quantities can be gauge 
dependent, however, all observable quantities should be 
gauge invariant. However this is not the case in the cold 
atoms experiment that the observable physical quanti- 
ties, e.g., the momentum distribution of the wave func- 
tion are gauge dependent. There is no paradox and on 
the other hand they are compatible in the case that the 
imaging of density of the expanding condensates in cold 
atoms experiment are in fact the canonical momentum 
of the original model [17j. This just reflects nothing but 
the exitance of vector gauge potential. 

In the case of synthetic magnetic field we choose the 
Landau gauge, which is : 

A = (0,27r/aO. (4) 

Obviously the hamiltonian in Eq. (j3]) is periodic in / with 
the period 1, thus we only need to study the properties 
of the system in the interval / G [0,1]. Here / = 
corresponds to the unfrustrated case and / = | is the 
fully frustrated condition. 

The wave function of the system in the lattice space 
ipi can be calculated using the Monte Carlo method. At 
low temperature the system undergo a superfluid transi- 
tion and there will be a peak in the central of the wave 



function in the momentum space, which is fourier trans- 
formed as: 

s 3 

where N s is the number of sites of the square lattice. 
The central peak of the momentum space describes the 
coherence of the phase of the condensates. It is analogous 
to the magnetization of a spin system, i.e., 

M=(ip ). (6) 

In cold atoms experiments, the observation of momen- 
tum distribution can be detected by the sudden release 
of the optical lattice. The absorption imaging is taken 
after a TOF period t. The density profile of the image 
can be written as : 

n(x) = (M/^) 3 |cD(k)| 2 G(k). (7) 

Here momentum k is related to position x by k = Mx/M 
under the assumption of ballistic expansion. u)(k) is the 
Fourier transform of the Wannier function and G(k) is 
the momentum space density matrix and is defined by: 

G ( k ) = ^E eJk ' (r, " rj) ^)- ( 8 ) 

S i,3 

For density coherent states at low temperature, the 
system forms a quasi-condensate where on each site of 
the lattice there is a small condensate and we can have 
(tpl^j) ~ iptipj- Thus the momentum space density ma- 
trix G(k) can be written as: 

G ( k ) «^E eik ' (r, " r3) ^' 

S i,3 

= N s -$ k \ 2 . (9) 

In this case, the momentum space wave function is 
related to the density matrix and then can be observed 
experimentally using the TOF imaging. In the following 
we will do the Monte Carlo simulation and investigate the 
transition behavior under different frustration /, which 
can be directed realized in cold atoms experiments. 

III. RESULTS AND DISCUSSIONS 

We now use the standard Metropolis Monte Carlo 
method with periodical boundary conditions to simulate 
the frustrated XY model. It has been used to investigate 
the BKT transition in Josephson arrays by Titel. This 
method is proved to always get believable results. The 
lattice sites are chosen as: L x L = 40 x 40. For each tem- 
perature and frustration of the system we use 10 7 Monte 
Carlo steps. 

Fig. [1] illustrates the evolution of expansion image from 
phase coherent ground state to the non-coherent state. 
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FIG. 1: An illustration of expansion image of the system for 
different temperature T at frustration / = 1/2. The BKT 
transition is taken place at about T — 0.5J/£;b as in fig.b. 
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FIG. 2: The central peak Go(G(k x = 0,k y = 0)) of expansion 
image as a function of temperature T for four different frac- 
tional frustration /. The insert is the corresponding image 
close to the ground state. The square points are the numeri- 
cal results with the error bar the standard deviation and the 
circle in each figure guides the critical transition temperature. 



At sufficient low temperature the system is in a super- 
fluid state and the density profile reveals regular sharp 
peaks at the momentum space lattices. That is due to 
the pairing of vortices. As the temperature increasing 
the peaks begin to decay and at the transition point at 
about Tc « 0.5J/kB the sharp peaks go to zero. That's 
the signature of BKT transition of 2D Bose gases. We 
can see that above the critical temperature there is disor- 
dered in the density profile of the momentum space and 
the system loses its coherent phase. 

In Fig. [2] we show the central peak Go(G(k x = 0, k y = 
0)) of expansion image as a function of temperature T 
for four different fractional frustration /. At T = the 
peaks have maximal values where the system is in a su- 
perfluid state with paired vortices and nonzero M. As 
the temperature is growing, the peaks are reduced and 
the vortices begin to unpair. The central peak drops to 
zero at the critical point and the transition temperature 
is guided by the circle. The inserts illustrate the ground 
state image for different frustration. For the fully frus- 
trated case in Fig. [2^ we can see the curve reveals differ- 
ent tendency as in Fig. [2)3 and Fig.[2t. The expansion im- 
ages in the latter case are the same as the non-frustration 
system. We note that ground state has 1/ f degenerate 
states while this can not be directly revealed from the 
expansion image. However, they can be recognized using 
printing phase technology [17j. 

We show the critical temperature T c at different frac- 
tional frustration / in Fig. [3j It reveals a unanalytical 
behavior of the transition temperature. We can see that 
there are peaks in some fractional magnetic field. The 



largest peak is the fully frustrated case in / = \ while 
the second largest is in / = |, which agrees with former 
experimental and theoretical work on superconducting 
Josephson arrays [7]. We note that the structure of the 
diagram depends on the geometry of the lattice. For 
a triangular structure the experimental and theoretical 
works reveal a secondpeak for the transition tempera- 
ture in / = j instead d, Q . 
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FIG. 3: Transition temperature T c of 2D BKT transition for 
different fractional frustration /. The critical point is corre- 
sponding to the absence of the central as illustrated in Fig. [21 
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IV. CONCLUSIONS 

In conclusion, we have studied the BKT transition of 
2D Bose gases in the synthetic magnetic field using the 
standard Metropolis Monte Carlo method. The criti- 
cal transition temperature is decided by the absence of 
central peak of the density matrix in momentum space, 
which can be detected in cold atoms by TOF expansion 
imaging. We have got the overall phase diagram of BKT 
transition for different fractional frustration. The ob- 
tained results agree with former experimental and the- 
oretical work on superconducting Josephson arrays (zj. 
Our results suggest that cold atoms can be used as an 
ideal system to explore the physics of the frustrated 2D 
systems. 

Note added: After this paper was almost finished, a 
preprint by Y. Nakano, K. Kasamatsu, and T. Matsui 



appeared on the arXiv.org [18]. These authors have stud- 
ied the finite-temperature phase structures of hard-core 
bosons in a two-dimensional optical lattice subject to 
an effective magnetic field based on the extensive Monte 
Carlo simulations. 
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